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QCD possesses a compact gauge group, and this implies a non-trivial topological 
structure of the vacuum. In this contribution to the Gribov-85 Memorial volume, 
we first discuss the origin of Gribov copies and their interpretation in terms of 
fluctuating topology in the QCD vacuum. We then describe the recent work with 
E. Levin that links the confinement of gluons and color screening to the fluctuating 
topology, and discuss implications for spin physics, high energy scattering, and the 
physics of quark-gluon plasma. 
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1. The Gribov copies and knots in the QCD ether 

In 1977, Gribov pointed out ^ the fundamental ambiguity in non-Abelian 
gauge theory due to the existence of multiple solutions of the gauge fixing 
condition, known today as “Gribov copies”. Much of the analysis was 
carried out in Coulomb gauge, but shortly afterwards Singer ^ showed that 
the problem persists in all other gauges, and thus no gauge fixing is possible 
in non-Abelian theories. Gribov has proposed to avoid the problem of copies 
by restricting the range of functional integration to small amplitude fields 
for which the solution to the gauge fixing condition is unique. He has 
demonstrated that such a restriction has a dramatic effect on the Coulomb 
interaction in QCD - it becomes confining, with Coulomb potential rising 
linearly with the distance, see for reviews. 

In this contribution (based on a recent work with E. Levin ®), we will 
take a different route - namely, we will argue that Gribov copies naturally 


1 
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arise from the compactness of the gauge group, and instead of eliminating 
them, we have to build a picture of the vacuum that accommodates their 
existence. The physical quantities in non-Abelian theories should thus be 
evaluated by summing over Gribov copies. When applied to the gluon 
propagator in QCD, this recipe leads to the confinement of gluons at large 
distances 

In fact, shortly after Gribov pointed out the existence of gauge copies, 
Jackiw, Muzinich and Rebbi ® showed that the copies, and associated with 
them discontinuities of gauge potentials are necessary to accommodate the 
vacuum configurations with arbitrary Chern-Pontryagin index. In an even 
number of space-time dimensions, such as our (3-1-1) dimensional Minkowski 
space, the Ghern-Pontryagin index is associated with topology-changing 
transitions. The discussion of topology of gauge theory is usually restricted 
to semi-classical topological solutions (instantons, monopoles, dyons, ...) 
and is thus in the domain of “non-perturbative QCD” traditionally ignored 
by the practitioners of perturbation theory. However a meaningful pertur- 
bartive expansion is based on fixing the gauge, and the existence of Gribov 
copies thus implies the impossibility to isolate perturbation theory from the 
topological nature of QCD. 

The relevance of topology in gauge theories stems from the compactness 
of the gauge group. Non-Abelian gauge theories always possess a compact 
gauge group, but even in the Abelian case we may choose to treat the U{1) 
group as compact - and in this case the theory acquires a non-trivial topo¬ 
logical structure, and the confined phase, see e.gJ. This fact suggests a 
link between topology and confinement. Compact gauge groups in general 
possess a homotopy map to the space-time manifold. For example, the 
map from the SU(2) C SU(N) subgroup of the gauge group to the Eu¬ 
clidean space-time sphere describes the classical instanton solution ®. In 
Minkowski space-time, instanton solutions represent the tunneling events 
connecting the degenerate vacuum states with different Chern-Simons num¬ 
bers 

X{t)= J (fxKo{x,t), (1) 

where Kq is the temporal component of topological current 

; (2) 

we have absorbed the coupling constant in the definition of the field. 
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The density of Chern-Pontryagin number mentioned above is the diver¬ 
gence of topological current (2) 

Qix) = d>^K^, (3) 

and is thus given by 

Qi^) = ( 4 ) 

The Chern-Pontryagin index is J d^x Q{x) = v\ for finite action field con¬ 
figurations V is an integer. As shown in Refs the presence of Gribov 
copies is necessary to accommodate configurations with an arbitrary (inte¬ 
ger) value of Chern-Pontryagin index. 

Prior to proceeding further, let us develop an intuitive picture of the 
vacuum structure implied by the topology of the compact gauge group. For 
simplicity, let us begin with the case of Abelian gauge theory - a familiar 
electrodynamics. In this case we limit ourselves to the first term in (2). The 
corresponding Abelian Chern-Simons number (1) has been known since 
Gauss’s work in the XIX century - it is a so-called “magnetic helicity” 
introduced in magnetohydrodynamics by Woltjer® and Moffatt^®: 

%{t) = J d^x A-B; (5) 

to be precise, the Chern-Simons number X{t) (1) is equal to 'H(t) divided 
by 167r^. 

Magnetic helicity measures the “knottedness” of the lines of magnetic 
flux. To quantify this statement and to characterize the topology of mag¬ 
netic flux, let us split magnetic field into N magnetic flux tubes with mag¬ 
netic fluxes (j)i. The magnetic helicity can then be written down 10 . 11 , 12,13 
in terms of the geometric invariant of the links between these flux tubes: 

N 

i,j 

where Si is the Calugareanu-White self-linking number, and Cij is the 
Gauss linking number®'. 


^The linking numbers in (6) do not always detect the topology of the link; the devel¬ 
opment of the appropriate universal knot invariants is a very active area in modern 
mathematics. 
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In infinite space, the solution of the Maxwell equations characterized by 
a non-zero value of (6) is very simple - it is a circularly polarized plane wave. 
The magnetic helicity is a parity-odd quantity, and the parity transforma¬ 
tion turns a left-handed circularly polarized wave into the right-handed one. 
In the presence of the boundaries, the solution^"^ becomes a Hopfion with 
linked lines of fields characterized by non-zero value of the second term on 
the r.h.s. of (6). 

The notion of magnetic helicity is extensively used in magnetohydrody¬ 
namics, where it is conserved in the absence of dissipation - in geometri¬ 
cal picture implied by (6), this means that the lines of magnetic field are 
“frozen” into the fluid and do not reconnect. Because the time derivative 
of magnetic helicity is 

^n{t) = ^ ( d^x A - B = -2 ( cPx E ■ B (7) 

dt dt J J 

from the definition of electric field E, we see that conservation of magnetic 
helicity implies J d?x E B = 0. In electrodynamics, this quantity is nothing 
but the Chern-Pontryagin index v defined above, since = —AE ■ B. 

The gauge field configurations with non-zero Chern-Pontryagin numbers 
(such as the ones signaled by the presence of Gribov copies in non-Abelian 
theory) thus have an important role: they induce reconnections of magnetic 
flux. 

In classical magnetohydrodynamics, reconnections of magnetic flux are 
often dramatic events - for example, in the corona of the Sun they are 
responsible for solar flares and coronal mass ejections spanning millions of 
miles (for review, see^®). In quantum theory, the physics during a magnetic 
reconnection can be described by adding the Chern-Simons term to 

the action of electrodynamics^®, where /is is the chiral chemical potential 
- since magnetic helicity AT® is parity-odd, /is is parity-odd as well. The 
equations of motion of the resulting Maxwell-Chern-Simons theory then 
exhibit the chiral magnetic effect (see for review and references) - 
the electric current J = e^/27r^ /is-B- The spectrum of photons in this 
theory differs from the one in usual Maxwell electrodynamics and exhibits 
a “chiral magnetic instability” - the photon energy acquires an imaginary 
part at soft momenta k < (e^/27r)/is reflecting a growth or decay of the 
photon fieldi®’ 20 . 2 i. 22 _ 

At high momenta k (e^/27r)/xs the photons behave as a free quanta 
of perturbation theory, but in the infrared region the behavior is totally 
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different - the photons either cease to propagate (are confined), or are pro¬ 
duced (this is a quantum analog of the solar flare induced by a magnetic 
reconnection). In the case of the vacuum, which is by definition the state 
with the lowest energy, the radiation of photons is clearly impossible and 
the physical choice for the imaginary part of the photon energy corresponds 
to the photon decay. If the magnetic reconnections take place in the vac¬ 
uum with a certain rate per unit time per unit volume, then this rate will 
determine the characteristic length over which the soft photon is allowed to 
propagate. We will see in the next section that this may be exactly what 
happens in QCD to soft gluons. 

The emerging picture of the vacuum filled by the knots is remarkably 
close to the one proposed by Lord Kelvin in 1867^^. Inspired by Helmholtz’s 
work on vortex motion in perfect liquid, Kelvin proposed that the observed 
atoms were in fact the knots of vortices in the ether. He wrote^^: ...the 
endless variety of [knotted or knitted vortex atoms] is infinitely more than 
suffeient to explain the varieties and allotropies of known simple bodies and 
their mutual affinities. While we now know that atoms are definitely not 
knots in the ether, Kelvin’s prophecy may still apply to the structure of the 
QCD vacuum. 

2. Veneziano’s ghost and fluctuating topology in the QCD 

vacuum 

The relevance of non-trivial topology in the QCD vacuum is revealed in the 
physical spectrum of the theory through the chiral anomaly^^’^® and the 
resulting non-conservation of the axial current: 

= 2NfQ{x) + ^(2tTO/)g/75g/, (8) 

/ 

where m/ are the masses of quarks, Nf is the number of flavors, and Q{x) is 
the density of Chern-Pontryagin number given by (3). The vacuum fluctu¬ 
ations in topology and the chiral anomaly thus lead to the explicit breaking 
of C/a(1 ) symmetry, and are responsible^®’^® for the anomalously large mass 
of the would-be Goldstone boson f, as we discuss below. 

Veneziano has demonstrated that the fluctuations of topology in the 
vacuum can be captured by introducing a massless “ghost” in the correla¬ 
tion function of the gauge-dependent topological current (2): 

K^^iq) =t j d^x {T{K^{x)K,{Q)}) ( 9 ) 
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where /i^ = Xtop is the topological susceptibility of pure Yang-Mills theory. 
Note that the r.h.s. of Eq. (9) has the “wrong” sign, which means that 
the ghost does not describe a propagating degree of freedom. Therefore 
the ghost cannot be produced in a physical process; however the couplings 
of the ghost (that describe the effect of topological fluctuations) can affect 
physical amplitudes. A similar “dipole” ghost had been earlier introduced 
by Kogut and Susskind in the analysis of axial anomaly in the Schwinger 
model. This procedure has been demonstrated to solve the Ua{^) problem 
in QCD 28.29 


The physical meaning of Eq. (9) becomes apparent if one compares it 
to the correlation function of the electron’s coordinate x{t) in a crystal, 
following Dyakonov and Eides 


i 


dt (T{a;(t)a;(0)}) 


u ;—>-0 
-)• 


1 

u;2 m* 


J_ d^E{k) 
a;2 dk"^ 


k=0 


( 10 ) 


where E = k"^/2m* is the energy of an electron with an effective mass m* 
and quasi-momentum /c in a crystal. The emergence of the pole in Eq. (10) 
signals the possibility of electron’s propagation in the periodic potential of 
the crystal by means of tunneling. It is important to note that the pole 
emerges not just from a single tunneling event (corresponding to the isolated 
instanton in QCD), but sums up the effect of many tunnelings throughout 
the crystalline lattice. 


The analogy between Eq. (9) and Eq. (10) becomes even more apparent 
if we choose the frame with q^^ = (w, 0) and use the analog of coordinate 
given by Eq. (1). The expression Eq. (9) then takes the form that is com¬ 
pletely analogous to Eq. (10): 


i 


dt 6*“* {T{X{t)X{0)}) 


d'^E{e) 

^ w2 w2 QQ2 


( 11 ) 


where V is the volume of the system, and E(0) = e(0)V is the energy 
of the Yang-Mills vacuum. The energy density of the vacuum e(0) is a 
periodic function of the 0 angle that is analogous to the quasi-momentum 
k in Eq. (10). At small 6, we can expand e{9): 

e(0)=8/r4^, (12) 


which exhibits the physical meaning of the quantity on the r.h.s. of (11 
- it is the topological susceptibility Xtop = of the Yang-Mills theory. 
Note that a term linear in 6 is forbidden by P and CP invariances of QCD. 
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The expression (11) clearly exhibits the role of the ghost - it describes 
the effect of topology-changing transitions characterized by non-vanishing 
Chern-Pontryagin indices in the vacuum. In the Abelian case presented 
in the previous section, the ghost would describe magnetic reconnections 
occurring with the rate (per unit volume, unit time) proportional to 

This interpretation becomes even more apparent if we consider the corre¬ 
lation function of the Chern-Pontryagin number density Q{x) (3) implied by 
(9), as we will now explain. The assumption of the ghost existence (9) im¬ 
plies the following Minkowski-space correlation function for Q{x) = d^K^-. 

i lim f d'^x (TlQ (x) Q (0)}) = < 0. (13) 

<1^0 J 

The quantity Q{x) is gauge-invariant, and excites from the vacuum physical 
states, the lightest of which is the flavor-singlet r/' meson'^. By writing down 
the spectral representation for the correlator (13) it is easy to see that the 
contribution of these physical states to the r.h.s. should be positive: 

z lim I d^x {T{Q (x) Q (0)}) = / ^ > 0. (14) 

How can we reconcile (13) and (14)? The key to the answer stems from 
the observation that in the physical world (in which the dispersion integral 
(14) is written), the theory possesses light quarks of mass m/, and in the 
chiral limit of massless quarks the correlation function of topological charge 
in (14) must vanish. This is because at m/ = 0 the topological charge is a 
full divergence of the gauge-invariant axial current (see (8)). 

For non-flavor-singlet axial currents that do not mix with gluon 

states, the vanishing of the correlation function of 9^ ^5 at zero momen¬ 

tum is assured by the Goldstone theorem, since the matrix elements of the 
divergence of axial current between the vacuum and the Goldstone boson 
(e.g. a pion) 

{0\df,J'^Ns\^)=^lf^ (15) 

vanish in the chiral limit of massless quarks as m,r 0. 

In the case of the flavor singlet axial current producing from the vac¬ 
uum the 7]' meson with matrix elements (0|J^|V) = \/‘2Nf frj' and 


'’For simplicity we do not discuss here the rj' — rj mixing. 
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ATsk) = \/2iV/ TO^/,r (normalization is chosen in accord with (8)), 
the situation is different because of the presence of topological fluctuations 
expressed by the relation (13). To construct the physical correlation func¬ 
tion in QCD, we thus have to add the r.h.s. of (13) and (14). Saturating 
the dispersion integral by the rj' peak and choosing in accord with (8) 
{0\Q\r]') = l/y/2Nf and ^(M^) = 4iV/ 2tz5{M'^ - m^,), we 

get 

* j ^ (0)})qcd = = 0: (16) 


from which we immediately get the Witten-Veneziano relation^®’^® for the 
rj' mass: 


m„ 


/2 = iNf = ANf 


d^e{6) 


96*2 


YM 


9=0 


(17) 


we have indicated explicitly that the topological susceptibility has to be 
computed for Yang-Mills theory (without light quarks). 


Subtracting a constant from the r.h.s. of (14) may seem strange - how¬ 
ever, we have to keep in mind that the dispersion integral in (14) diverges 
in the high mass region that we have omitted, since on dimensional grounds 
p(M2) ~ (as reproduced by the gluon loop) - therefore, it requires sub¬ 
tractions, and this is what the term represents'^. The Witten-Veneziano 
prescription described above thus implies that the correlation function (13) 
should contain a contact negative term, (rjQ (x) Q (0)}) ~ —/i^ (5‘*(x). 
There is an ample evidence for this contact term from the lattice simula¬ 
tions, see Fig.l from (the correlator shown in Fig.l is Euclidean, hence 
the sign difference). 


The physical interpretation of the contact term is suggested by our 
discussion in section 1 - it describes the “white noise” of topological re¬ 
connections, with strength given by the topological susceptibility Just 
like the introduction of white noise in kinetic theory allows one to describe 
the effect of the medium on a test particle, the ghost allows to capture the 
impact of fluctuating topology. Compared to “microscopic” approaches to 
QCD vacuum, such as the instanton liquid model and its descendants (for 
review, see ^^), the approach described here does not identify the dynamics 
of topological transitions - instead, we parameterize the rate of transitions 


The T-ordered product in (14) has to be understood as a Wick, and not Dyson, one®®. 
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by a single dimensionful parameter, topological susceptibility. On the other 
hand, it does not rely on a (questionable) semiclassical approximation. 

Usually, the discussion of topological effects is limited to the U^(l) 
problem and correlation functions of color-singlet correlators, as presented 
above. However, the fluctuating topology should also affect the propagation 
of gluons, similarly to the case of photons discussed in the previous section. 
We will now see that this is indeed the case. The Veneziano ghost defined 
by (9) will allow us to address this problem in QCD without resorting to 
the complicated microscopic dynamics of topological transitions. 


Normalized tofmlogical charge density correlator 



Figure 1. The lattice QCD data on the correlation function of Chern-Pontryagin charge 
density at three lattice spacings a; from^^. 


The observation made in our recent paper^ with E. Levin is that (9) and 
(2) define an effective ghost-gluon-gluon vertex T^ {q,p)- Since the topolog¬ 
ical current (2) is a pseudo-vector, the ghost is a pseudo-vector “particle” - 
so the interaction of a gluon with, say, a left circular polarization with the 
ghost will transform it into a right circularly polarized gluon. Because the 
QCD vacuum is P and CP even, the ghost can only appear in the loops. 
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Using this vertex, we can re-write the correlator Eq. (9) at small as 
follows (see Fig. 2-a): 

1 /* 1 ^ 

J ( 18 ) 

from which we find that 

^4 

iq,p) iq,p) oc , for g < p. (19) 

The vertices T^ {q,p) describe the interactions of soft gluons with topolog¬ 
ical fluctuations. 

Indeed, the gluon propagator is now the solution of the Dyson-Schwinger 
equation with S (p) given by 


/a {q,p) (foP) (20) 

where is the propagator of the ghost. In evaluating the integral of 
Eq. (20) we assume that q p since E^’s describe the non-perturbative 
effects at small momenta. Therefore, 

dir 1 

^iJ.y{p) - TT-m / dq^q^ T^j,{q,p)^^T^{q,p) = (21) 

(27r)4 q^p^ 

where we used Eq. (19). 


We can now write down the Schwinger-Dyson equation for the gluon 
propagator (p) = G (p) : 

Gip) = \ + \j:{p)G{p) (22) 

pZ pZ 


with the solution 


G{p) 


1 

p2 - S (p) 


1 



(23) 


The propagator (23) has some remarkable properties. First, G {p) has 
no gluon pole in the physical region. Instead, this propagator has only 
complex poles at = zbz/x^. As a result, gluons cannot propagate at large 


^This is the expression in Feynman gauge; since the self-energy does not change 
the polarization of the gluon, the Dyson-Schwinger series can be re-summed in the usual 
way. 
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distances - in other words, they are confined. Second, the propagator of 
the type of Eq. (23) as shown by Gribov ^ and Zwanziger eliminates the 
gauge copies. Note however that we have obtained this propagator not by 
restricting the range of functional integration to the region where the copies 
do not appear, but by taking account of all copies into account through the 
topological fluctuations. 

Let us illustrate the meaning of our result by appealing again to the 
analogy with the crystal, see (10). The Gribov copies in this example are 
analogous to the electron states localized around different sites in the crys¬ 
talline lattice. When we discuss the wave function of the electron around 
one of the sites, in general we cannot ignore the presence of the tunnel¬ 
ing throughout the lattice. The tunneling can be ignored only if we limit 
ourselves to the tightly bound electron levels. However the tunneling can 
be consistently taken into account if instead of an individual bound state 
we consider the electron bands. Electron energies in these bands acquire a 
finite width, similarly to what we have observed for the gluons. 

Our propagator Eq. (23) results from the admixture of the ghost to the 
perturbative gluon (see Fig. 2-b), with an amplitude defined by the Yang- 
Mills topological susceptibility In paper^, we proposed to call this 
coherent mixture of a gluon and a ghost a “glost”. Unlike the ghost, the 
“glost” can be produced in a physical process, but unlike the perturbative 
gluon, it is confined and cannot propagate at distances larger than ~ ~ 

1 fm. 



Figure 2. Fig. 2-a : Eq. (9) in momentum representation; helix lines represent gluons 
and the dashed line depicts the ghost. Fig. 2-b : the gluon dressed by the interactions 
with the ghost: a “glost”, see Eq. (20) for the corresponding self-energy expression. 


At distances ~ ~ 1 fm, the use of the glost propagator also gives 

rise to the screening of color charge^, leading to the freezing of the effective 
coupling in the infrared IR limit for pure gauge theory, or to the vanishing 
of the coupling in the IR in QCD with light quarks, see Fig. 3. It has been 
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shown by Dokshitzer that the experimental data indicate the relative 
smallness of the QCD coupling in the IR region: 

ao = — / dkas{k) « 0.5 for fij = 2 GeV. (24) 

k'l J 

In our approach we get Oq — 0.59 for renormalization point k = Mz, in 
reasonable agreement with (24). 





Figure 3. The running constant 05 as a function of momentum k in our model for QCD 
with light quarks (red line that goes to zero at k —>■ 0) and gluodynamics (blue line); 
from®. 


Note that close to the deconfinement transition, the topological suscep¬ 
tibility vanishes reflecting the restoration of C/zi(l) symmetry 34,35^ Since at 
p, —>■ 0 the gluon propagator becomes perturbative, the restoration of 17^(1) 
symmetry and deconfinement should occur at the same temperature in 
agreement with the lattice data however close to Tc the non-perturbative 
interactions induced by topological fluctuations are important. It would be 
interesting to explore the implications of this approach to QCD thermody¬ 
namics, and to the rates of electromagnetic emission from the quark-gluon 
plasma - this work has already begun 

It is of fundamental interest to establish the microscopic dynamics re¬ 
sponsible for the fluctuating topology captured by the ghost. A recent study 
within the “deformed QCD” attributes these fluctuations to the topological 
order in the vacuum The imaginary part of the gluon self-energy signals 
an instability w.r.t. decay to the true ground state - what is the structure of 
this ground state? The clue to the answer may be provided by the Abelian 
example considered in section 1 - there, the true ground state of the photon 
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field is the (self-)linked knot of magnetic flux. The knots repeatedly emerge 
in studies of the ground state of non-Abelian gauge theories'*^’^^’^^, but it 
is clear that much remains to be done. 

In the meantime, the simplistic approach described above captures the 
link between the confinement and the fluctuating topology in the QCD 
vacuum, and provides a practical way of computing power-suppressed cor¬ 
rections to QCD amplitudes. In particular, the glost propagator leads to 
the exponential fall-off of the high-energy hadron scattering amplitude at 
large impact parameters needed to satisfy the Froissart bound; this can 
solve the long-standing problem of the perturbative approach in describing 
high energy scattering In QCD amplitudes the coupling to the ghost 
flips the gluon and (through loops) quark polarizations, and can give rise 
to spin asymmetries that are different from the usual perturbative ap¬ 
proach. 
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